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Abstract. Gaussian stochastic diffusion processes are used to derive cosmic mass functions. To get analytic rela- 
tions previous studies exploited the sharp fc-space filter assumption yielding zero drift terms in the corresponding 
Fokker-Planck (Kolmogorov's forward) equation and thus simplifying analytic treatments significantly (excursion set 
formalism). In the present paper methods are described to derive for given diffusion processes and Gaussian random 
fields the corresponding mass and filter functions by solving the Kolmogorov's forward and backward equations in- 
cluding nonzero drift terms. This formalism can also be used in cases with non-sharp fc-space filters and for diffusion 
processes exhibiting correlations between different mass scales. 
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1. Introduction 

The Press-Schechter (1974) formula of the cosmic mass function has found many important applications to basic 
problems in modern cosmology (e.g., Narayan & White 1988, Cole & Kaiser 1988, 1989, Henry & Arnaud 1991, White 
& Frenk 1991, White et al. 1993, Eke et al. 1996, Mo et al. 1996, Mo & White 1996, Matarrese et al. 1997, Mathiesen 
& Evrard 1998, Borgani et al. 1999). Although the approach appears questionable when the theoretical masses are 
compared with the group finder masses in detailed N-body simulations (e.g., Bond et al. 1991, White 1996) the overall 
statistical agreement is found to be sufficiently good (Efstathiou et al. 1988, Lacey k, Cole 1994, but see also Sect.||). 
Improved versions relax the various assumptions of the original derivation (Lucchin & Matarrese 1988, Bond et al. 
1991, Lilje 1992, Cavaliere & Menci 1994, Jedamzik 1995, Monaco 1995, Kitayama & Suto 1996, Yano, Nagashima & 
Gouda 1996, Valageas k Schaeffer 1997, Chiu et al. 1998, Lee & Shandarin 1998). 

Apart from providing a convenient fitting equation, the excursion set formalism in the form presented in Bond et 
al. (1991) serves as the basis for several fundamental theoretical developments aiming for a full treatment of structure 
formation in a hierarchically clustering universe. Extended theories for hierarchical clustering (Bower 1991), biasing 
schemes (Mo & White 1996), and merging histories (Lacey & Cole 1993, Kauffmann & White 1993) are developed or 
can also be derived in the end by assuming the validity of the excursion set arguments introduced by Bond et al. to 
derive the Press-Schechter formula with the correct normalization and several related conditional probability densities. 

For the derivation of the Press-Schechter function the excursion set formalism assumes that at a fixed spatial 
location the change of the density contrast successively smoothed with a decreasing filter scale, R, can be described 
as a sample path of an abstract diffusion process (Bond et al. 1991, see also Peacock & Heavens 1990). The counting 
of subclumps (clouds-in-clouds) is avoided by computing the rates at which the sample paths meet at their largest 
mass extent a scale-independent critical density contrast barrier, / c . For this process Bond et al. derived a general 
Fokker-Planck (Kolmogorov's forward) equation consisting of a drift and of a diffusion term. The application of the 
sharp fc-space filter condition yields a zero drift term so that the resulting much simpler Fokker-Planck equation can 
be solved easily and leads to the Press-Schechter mass function with the correct normalization. 

Unfortunately, analytic results are obtained with the excursion set formalism only within the theoretical framework 
of Markovian processes (see below) and for the sharp fc-space filter with a mass assignment scheme which clearly 
contradicts numerical experiments on the halo-by-halo basis. It is the basic aim of the present investigation to extend 
the excursion set formalism to non-sharp fc-space filters but without leaving the framework of Markovian random 
processes and thus preserving the simplicity of the original approach. In the present context, Markov processes are 
assumed to be driven by deterministic and stochastic 'forces', and have the simplifying property that when the process 
states, in form of amplitudes of filtered density contrasts, are known on a given mass scale and larger, the process 
states on smaller scales depend only on the actual scale and are independent of the larger scales - it seems suprising 
that under this Markov condition, e.g., useful biasing schemes can be derived, but see below. 

The present paper discusses both the Kolmogorov's forward and backward equation including nonvanishing drift 
terms. Analytic solutions are presented for several interesting cases offering the possibility to derive cosmic mass 
functions independent of the sharp fc-space filter assumption. It will be shown that certain Gaussian stochastic diffusion 
processes lead to more realistic filter profiles and to correlations between different mass resolution scales which are 
expected to have important implications on various progenitor statistics and biasing schemes. 
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Sectiong gives an informal presentation of the basic concepts of the excursion set formalism. Section^ summarizes 
the formal properties of the relevant diffusion processes and their sample paths. The connection between diffusion 
processes and mass functions is outlined in Sect.|j. The relation between the basic pseudo time (resolution) variable 
and the filter radius is discussed in Sect.||. In Sects. ^| and [7] the methods are applied to order-0 (Wiener) diffusion 
processes, characterized by drift coefficients which are independent of mass and filtered density contrast, and to order- 1 
(Ornstein-Uhlenbcck) diffusion processes, characterized by drift coefficients which are also independent of mass but 
which scale linearly with filtered density contrast. For both processes the corresponding mass and filter functions are 
derived. The results are summarized and discussed in Sect.|^. Appendix A gives an overview of the basic properties of 
order- 1 diffusion processes. 

Forthcoming papers will analyse the new mass assignment scheme, and will compare the mass functions presented 
here and others deduced from more complicated stochastic diffusion processes with detailed N-body simulations and 
with observed mass and luminosity functions of rich clusters of galaxies obtained with, e.g., the REFLEX cluster 
sample (Bohringer et al. 1998, Guzzo et al. 1999), and the HIFLUGCS cluster sample (T.H. Reiprich & H. Bohringer, 
in preparation). 



2. Excursion set formalism 

The differential mass function, n(M), gives the number of virialized objects of mass M, per unit mass and (comoving) 
volume element. We want to derive n(M) using the field of density contrasts, f(r), at the comoving coordinates, r. 
Here, / is defined by the ratio f(r) — p(r)/p— 1, where p(r) is the comoving mass density at r, and p the mean 
comoving mass density. The following notation is simplified by neglecting any evolution of / with cosmic time. However 
the final equations can be transformed to any epoch by multiplication with linear growth factors. For the derivation 
of n(M) it is assumed that the location where an virialized object of mass M will be formed corresponds to the 
location where /(r), filtered with the radius R(M) to give f(r, R), has reached a universal critical threshold value, f c . 
The density field is divided into infinitesimal small volume elements. A sequence of filtered three-dimensional density 
fields is constructed with gradually decreasing filter radii. The excursion set formalism identifies the cumulative mass 
fraction cx n(> M) with the fraction of volume elements in which the density contrast lies above f c when smoothed 
with any filter of radius greater than or equal R(M). 

At a fixed spatial location, r, the sequence of filtered density contrasts, {f(r, R)\R — 00 . . . 0}, can be regarded as a 
one-dimensional trajectory. The statistical description of an ensemble of such sample paths gets especially simple when 
they result from a Gaussian random process with independent increments, i.e., from a Wiener process. For filtered 
density contrasts this is the case when f(r) is a Gaussian random field uniquely described by the fluctuation power 
spectrum, P(k), as a function of comoving wavenumber fc, and when the top-hat filter in Fourier space, i.e., the sharp 
fc-space filter (see Eq. || below) is used: 

The power spectrum of the density field is defined via the Fourier transform of the unfiltered density contrasts, 
/(fc) = FT[/(r)], by the expectation value E[/(fc) /*(&')] = ( 27r ) 3 P(fc)<fo( fe ~ where S v is the Dirac delta dis- 
tribution and * assigns the complex conjugate. In this case the discrete change in f(r, R) caused by the transition 
fc/v, — ► fc_R + Afc# of the maximum comoving wavenumber, kn, which passes the filter is a Gaussian random variable 
(each change of f(r, R) when R is changed comes from the inclusion/exclusion of a new independent fc shell) with the 
variance Act 2 = <r 2 [l/(fc/f + Afc^)] — a 2 [l/kn] 1 where 

° 2 ( R ) = 7o~TT / d 3 k\W R (k)\ 2 P(k). (1) 



(2tt) 3 Jo 
The sharp fc-space filter, 

and Eq. (Q) define (for spherically symmetric filters) a general resolution variable, 

A(R) = — / dkk 2 P(k). (3) 
Jo 

For this choice, the /(A) values perform some kind of Brownian motion with A(i?) as a pseudo time variable. Notice 
that for hierarchical clustering, A as well as a 2 increase with decreasing mass and filter scale. For Afc^ — > a trajectory 
corresponds to one realization of a continuous limit of a simple random walk: the Wiener process. 

Hence, the behaviour of these trajectories are governed by a simplified Fokker-Planck diffusion equation which 
gives the probability density, n, to find trajectories with f(r, R) at A. In this picture the differential mass function 
can be obtained from the rate at which the random trajectories meet for the first 'time' the absorbing barrier f c at a 
given mass or A[R(M)] scale. The final problem is to get solutions of the Fokker-Planck equation under the boundary 
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conditions that all sample paths start at A = with / = and that absorption of a sample path at A corresponds to 
n(/ = / c , A) = (see Sect. | for a more detailed discussion of these boundary conditions). 

The excursion set formalism in the form presented above provides a simple mean to deduce the Press-Schechter 
mass function. The major reason for the simplicity of the prescription results from the application of the sharp fc-space 
filter. The choice of this filter appears to be rather ad hoc and not thoroughly motivated. Its strongly oscillating shape 
in configuration space does not seem to be optimal to frame the region of primordial material that ultimately collapses 
to form a virialized halo. Moreover, the absence of any covariances between different mass scales is not expected for 
realistic physical systems. In the following a more general (though more formal) description of the evolution of II is 
given, by considering a general diffusion process which offers the possibility to study also non-sharp filter functions 
and the resulting mass functions. It will be seen that certain diffusion processes lead to nonzero covariances between 
different mass scales without leaving the general framework of Markovian random processes. These processes are thus 
expected to yield more realistic formation histories and biasing schemes. 



3. Gaussian A-/ diffusion processes 

A Gaussian random field of the mass density contrast, /, is filtered on a continuous set of decreasing scales R, giving 
at a fixed spatial location in configuration space a continuous sequence of filtered density contrasts, /, at the scale 
R, providing a measure of the spatial resolution, A(R). For scale-invariant power spectra, P(k) oc fc™, the resolution 
variable has the form A oc i?~("+ 3 ) (Eq. ||). The diffusion process is restricted to a continuous - A, continuous - / 
state space, where A and / assume the roles of the usual time and spatial coordinates, respectively. In general it is 
expected that the sample paths of the diffusion process start with / = /o = at A = 0. Formally, a Gaussian A-/ 
diffusion process corresponds to a Markov process for which the transition probabilities, n(/', A + AA, /, A), with the 
sharp value / at A and the random variable /' at A + AA, satisfy a continuity condition (similar to Eq. (a) in Arnold 
1974, p. 40), and 

^ m n AT / d ? (/ - I')™ A + AA ' A ) = A ™V-> A ) + °( £ ) W 

AA^O AA J\f'-f\< e 

exists for e > 0, m = 1,2, and A m — for m > 3. The A\ coefficient describes the deterministic drift or the 
instantaneous rate of change of the mean, and A2 the stochastic diffusion or the instantaneous rate of change of the 
squared fluctuations of the sample paths giving the expectation, E[/(A) — /(A')] m = A m (f, A)AA + O(AA). Note 
that the Markovian independence property introduces the simplifications mentioned in White (1997) concerning the 
application of the extended Press-Schechter theory to formation histories and biasing schemes: Because the Markov 
property is used to derive extended Press-Schechter, a say 10 12 M halo at redshift z = 1 does not 'know' whether it will 
be incorporated into a rich cluster or remain in a void at z = 0. The galaxy populations in protocluster and protovoid 
halos will thus be the same. However, the present investigation concentrates on a wider class of diffusion processes 
than studied by the extended Press-Schechter approach, some of them exhibiting nonzero correlations between different 
mass scales. The expectation E[/(A) — f(A')] m leads with the inclusion of some additional heuristic arguments (see 
the formal treatments in, e.g., Friedman 1975, Kloeden & Platen 1995) to the Ito stochastic differential equation, 



d/(A) = -A 1 (f,A)dA+y/A 2 (f,A)dw(A). (5) 
The second term on the right-hand side is the fundamental or Wiener stochastic process, i.e., a Gaussian process 
with the independent infinitesimal increments, dw(A) — w(A + dA) — w(A) = K(A)dA. A prescription for the 
construction of w(A) is given by eqs. (§) and @. The stochastic forces', K(A), are Gaussian with infinite variance 
and no correlations between different A shells. They drive the /-A diffusion as illustrated by the Langevin-type 



equation, df /dA = —A\(f,A) + y A-z(f, A) K(A). The increments, df ', are of the order of the standard deviations 

of the fundamental process, that is, df /dA V AA/AA, which tends to infinity as A A — > (see also Wiener 1930, 
Doob 1942). The Ito formalism is used to handle the resulting almost surely continuous but almost surely nowhere 
diffcrentiable sample paths in a mathematically well-defined way. An alternative would be the application of the 
Stratonovich formalism, but the results derived here do not depend on the specific choice. 

Closely related to the stochastic diffusion Eq. (|^) are Kolmogorov's forward (Fokker-Planck) equation, 

dU(f ,AoJ,A) = £( / ;A) n(/ , Ao ,/,A), £(/,A) = A) + \^Mf, A) , (6) 

with the forward variables, / and A, and Kolmogorov's backward equation, 

an(/ A ,/,A) = _ £ t ( / 0)Ao )n(/o,Ao,/,A), £t(/ ,A ) = -A 1 (/o,A )^ + iA 2 (/ ,A )^, (7) 
CA df 2 dffi 
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Fig. 1. Sketch of the large-scale behaviour of sample paths of low-order stochastic diffusion processes. Continuous lines: sample 
paths of a Wiener process with no drift term and the diffusion parameter A2 = 1. Dotted lines: sample paths of diffusion 
processes with the same diffusion parameter but nonzero drift terms; left panel: order-0 process with A\ — —0.5 (Eq. ^|). Right 
panel: order-1 process with A\ = 0.2/. 



with the backward variables, /o and Ao Here, & is defined as the formal adjoint of the elliptical operator C (e.g., 
Kloeden & Platen 1995, Sect. 2.4, Arnold 1974, Sect. 2.6). In Eq. (g) the backward variables / and A , and in (jf)) 
the forward variables / and A are essentially constant and enter only through boundary conditions. The sign of A\ is 
chosen to be consistent with (||) . The Ito formula (see Eq. [A.2j in Appendix A) can be used to verify the equivalence of 
the Fokker-Planck equation and (||). For /-A diffusion processes Bond et al. (1991) showed that Eq. (^) can be derived 
via the Chapman-Kolmogorov equation in the standard way. For these diffusion processes it is also straightforward 
to show the validity of Kolmogorov's backward equation. Hence, (|^) can be used also for more general /-A diffusion 
processes. 

The Ito stochastic differential equation cannot be solved in general. However, as a guideline consider equilibrium 
solutions with dH/dA = 0. Simple solutions of @ exist if A2 is independent of / and A, that is, A2 = D = const. For 
this choice the correct form of A\ is needed which is consistent with steady state. For A\ independent of A, A\ < 
would lead to / values 'accelerated' to large |/| with increasing A. Additional dependencies of Ai on even powers of 
/ would not change the situation. However, if Ai depends on odd powers of /, Ai > would lead for positive / to be 
less positive and for negative / to be less negative. Expanding A\ in odd powers of / thus yields stable deterministic 
states. In this respect important cases are the order-;/ processes, 

df(A) = -Af"dA + V15 dw(A) , (y = or odd) . (8) 
For v > 3, hitherto only small-noise expansions are known (e.g., Gardiner 1997, p. 185). We are thus left with the 
Ansatz v — and v = 1. Some historical notes on the resulting equations can be found in Wax (1954). Although it is 
not possible to give graphical presentations of sample paths of stochastic diffusion processes, Fig.|l| tries to compare 
trajectories of order-0 and order-1 diffusion processes obtained with 



/(A,) = -^[/(Ai-OJA, + ^/A 2 [f(A i . 1 )} Wl (i = l,...,N), (9) 
the ith resolution variable, Ai, and Wi from Wiener's construction as described in Ito & McKean (1974), 



2™-l 



'7T * ' * ' V 7T 

n=l k=2 n ' 



A ^o + E E f-^T^9 k (0<Ai<n), (10) 



where g^, k > are independent and Gaussian distributed with zero mean and variance 1. For m — > 00, Eq. ( J10|) 
converges uniformly in < Aj < tt with probability 1 and is a standard Brownian motion. The sample paths in Fig. Ill 
are computed with m = 10 and N = 1000 in the convergence range (larger m values do not significantly improve 
the figure). For the order-0 process we set A 1 [f(A i -i)] = —fi, for the order-1 process Ai[/(A,_i)] = Af(A i ^ 1 ), and 
A 2 [/(Ai_i)] = 1 for both (see Fig.|l| for more details). Fig. |l| illustrates the exceedingly irregular motion of order-0 
processes with an almost infinite 'speed', when / and A are regarded as the position and time of a Brownian particle, 
respectively. Order-0 processes have no stationary distributions, i.e., as A — > 00, any finite point moves to infinity with 
probability 1. Order-1 processes include a linear drift term, which acts like dynamical friction, and which yields the 
only Gaussian Markov process in one real variable with a stationary distribution. The process is ofte n use d to model a 
more realistic Brownian motion, in which A and A' are significantly correlated if |A — A'| ~ 1/A (Eq. A. 4 in Appendix 
A). 
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Fig. 2. Comparison of sample paths obtained for a high-density region (upper two trajectories), for a low-density region (lower 
two trajectories), and for an average-density region (middle two trajectories) by filtering a one-dimensional Gaussian random 
field (spectral index n — 0) using a sharp fc-space filter (continuous lines) and a truncated exponential filter (dotted lines) with 
A = 0.5. The filter scale, R, is set to arbitrary physical units. The systematic differences between the trajectories obtained with 
the sharp fc-space filter and with the truncated exponential filter reflect the differences between the trajectories of order-0 and 
order-1 diffusion processes shown in Fig.[j]. 



In Sect.|J it will be shown that each filter used to smooth Gaussian matter distributions and to extract masses of 
isolated objects is closely related to a specific diffusion process. To illustrate this point, sample paths are computed 
for specific locations (a high, average, low-density region) within a one-dimensional Gaussian random field which 
is smoothed with filter functions where the corresponding diffusion processes are known. Fig. || shows this for the 
sharp fc-space filter giving sample paths which correspond to order-0 diffusion processes (see Sect.[|), and for the 
truncated exponential filter giving sample paths which correspond to order-1 diffusion processes (see Sect.J^). Despite 
the artificial smoothing of trajectories at large 1/R caused by the limited numerical resolution in this scale range, the 
sample paths of the mathematical models of order-0 and order-1 diffusion processes shown in Fig.^ and the sample 
paths obtained from filtered Gaussian random fields obtained with the process-specific filter functions shown in Fig.|^ 
look qualitatively quite similar. It is thus seen that order-0 and order-1 diffusion processes are supposed to give a 
realistic description of the complex behaviour of the density contrasts at a given location in a Gaussian random field 
when filtered with process-specific filter functions and different filter widths. 

Concerning the physical meaning of the diffusion and drift parameter it is important to realize that whereas the 
diffusion parameter basically determines the relation between immediately adjacing mass scales, the drift parameter 
determines the relation between largely differing mass scales. 

Ignoring for the moment the deterministic drift subprocesses, the drift-less diffusion process is assumed to be 
independent of mass and density contrast. The resulting random behaviour of the filtered density contrasts is thus 
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determined by the inclusion or exclusion of a new independent k shell. Because of the specific choice of the resolution 
variable (see Eq. ||) and the normalization of the filter (see Sects. 6.2 and 7.2), the variance of the Gaussian noise 



process and thus of the value of the diffusion parameter, A2, is per construction always unity. With the addition of 
a deterministic drift subprocess the idealistic pure diffusion process can be converted into a more realistic process 
because the idealized assumptions of the Wiener process can be relaxed. 

The physical meaning of the drift parameter becomes clear from the deterministic motions induced (see Figs. [I] and 
^). As an example, for A-independent drift parameters (the diffusion processes discussed here have this property) the 
potential 

U(f) = [ A 1 (f')df' (11) 



might be defined (e.g., Risken 1984, Chap. 5), leading for the order-1 process with A > to a parabolic potential which 
restricts the otherwise freely streaming trajectories to smaller values of |/| (see right panel in Fig.|l|). In Appendix A 



(Eq. A. 4) it is shown that the corresponding 'external deterministic force' introduces realistic nonzero covariances 
between different /(A), that is, it restricts the trajectories to a smaller / interval around / = which has the following 
two relat ed c onsequences: (1) It results in a weakening of the sharp fc-space filter at larger comoving wavenumbers k 
(see Sect. |7.2| ) . This is of great cosmological significance because the corresponding filters in configuration space are 
comparatively smooth (see Fig. |j) and thus frame a more realistic region of material that immersed into a particular 
halo. (2) It results in a correlation between different mass scales and thus effectively weakens the strict Markovian 
property, however, without leaving the general framework of Markovian processes. These correlations are expected to 
give more realistic formation histories and biasing schemes compared to the original Press-Schechter approach. 

On the other side, order-0 processes with lead to systematic shifts of the trajectories in the direction of smaller / 
or larger values (left panel in Fig. Fy). This might be interpreted as a A (mass)-dependent change of the critical density 



threshold, / c , which will be discussed in more detail in Sect. 6.1 



4. Diffusion processes and mass functions 

In the following the general relation between diffusion processes and mass functions is developed. As already stated 
in Sect.^ the excursion set formalism assumes that the cumulative mass fraction n(> Af)/n to t, where n to t is some 
normalization constant, is identified with the fraction of volume elements above f c when smoothed with a filter of 
radius > R(M) which is directly related to A(R) (see below). At least for the high masses it seems reasonable to 
associate the change in total volume of the excursion region with the net mass in objects of this scale (see Bond et al. 
1991, and the useful graphical visualisations given therein). 

In the language of diffusion processes described in Sect.||, the differential mass function is proportional to the 
probability that a particular trajectory is absorbed at the threshold f c in a given A, i.e., mass range. It is thus 
necessary to investigate the first passage A distribution of the trajectories supposing an absorbing barrier at f c > fo- 
Let II(/o, /, A) be the probability density to have trajectories at /(A) and that no trajectory has reached f c within the 
interval (0, A). Notice the difference in the definition of II here and in Sect.||. Hence, the probability that absorption 
has not yet occured at A is (sometimes the fo variable will be omitted in the argument lists of F, II, and g) 

F(/o,/ c ,A) = t n(/ ,/,A)d/. (12) 



The lower limit, / = — 1, must be extended formally to — 00 in order to fulfill for all A or mass scales under consideration 
the assumption of a Gaussian diffusion process. Notice that F(fo, f, A) and U-(fo, f, A) satisfy (H@). The probability 
density function of the first passage of trajectories at A is determined by the 'loss-rate' of sample paths at A, i.e., 

5 (A,/ c ,/o) = - — y Yl(f J,A)df. (13) 

In the following g(A, f c , fo) is called the pseudo mass function. The probability density function g(A, f c ,fo) can 
be determined by solving Kolm ogo rov' s fo rward equation ([|) for II(/o,/, A) under appropriate boundary conditions 



(mirror image method, sec Sect. and 7.1). For A-independent f c , a second method determines the Laplace transform 
of g(A,f c Jo), 

S*(A,/c,/o) = / e- XA g(A,f c J )dA, (14) 
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and solves the Laplace-transformed Kolmogorov's backward equation (0) in conjunction with appropriate boundary 
conditions (Laplace transform method, see Sects. 6.1 and 7.1). The most difficult part of the latter method is the 
inversion of dl4). To obtain the differential mass function. 



n(M) = ^ 
M 



dA(M) 



dM 



<?(A,/ C) /o) 



(15) 



we have to derive A(M) which is closely related to the filter function of the diffusion process. For Gaussian matter 
distributions, this relation is obtained in the following way. In the first step, the equivalence between the mass variance 
(0) and the variance of the /'s, 

J d 3 k\W R (k)\ 2 P(k) = var[/(A)], (16) 

is used to derive the Fourier transform of the filter function, W R (k). The variances of the f's can be obtained with 
the Ito formalism in a straightforward way (see, e.g., Appendix A). In the second step, the resulting window function, 
W R (k), is transformed into configuration space to give W R (r), from which the filter volume is determined by V(R) = 
1/Wr(t = 0). The mass associated with the filter is then obtained from M(R) = pqV(R), where po is the present 
mean mass density of the Universe. Notice that this step follows the general conventions as introduced in Lacey & 
Cole (1994). In the last step, the general relation between A and R (see Eq. |l^) is applied to give the wanted A(M) 
relation. Sect.|5| presents a discussion of the general A(R) relation. 



5. The A-R resolution scale relation 

Within the framework of stochastic diffusion processes, mass functions are determined by the probabilities that the 
continuous sample paths of the diffusion process are first absorbed at / = f c . Analytical results are obtained by 



E[f(R).f(R')} = — 3 / d 3 kW R (k)W R ,(k)P(k) = <r 2 (max{R,R!}). (17) 



assuming no correlations between / and its infinitesimal increments, df, but (see, e.g., Eq. A. 8 for A = in Appendix 
A) 

1 

(2V) 3 

The last equality can be achieved if the filter function, W R (k), has the form of the sharp fc-space filter (Eq. ||), and 
the resolution variable, A, is defined by A = a 2 (R). This is the prescription of the fundamental stochastic diffusion 
process. Hence, the A-dependence of the probabilities II(/o,/,A) can be obtained from the simplified Fokker-Planck 
equation, 

A n(/ , /, A) = n(/ , /, A) , (is) 

with the diffusion coefficient, D = 1. It is a classical exercise of the Ito stochastic formalism to show that ( [T8| ) can be 
deduced from (Q) and vice versa for the constants A\ = and A% = D. Although these assumptions (including boundary 
conditions, for example, eqs.[24|) suffice to derive the exact form and normalization of the Press-Schechter function, 
they are for several reasons not very realistic. Zero-correlations between / and its increments are only attainable if 
the stochastic force, K, has zero mean and the covariance E[i£"(A') ■ AT (A)] = cov(A' — A) = D Sd(A' — A). This 
gives a white noise power spectrum and thus an infinite variance. As A — > 0, solutions of the simplified Fokker-Planck 
equation tend to the Dirac delta distribution. As A — > 00, the probabilities, II(-), tend to zero, and no stationary 
solution exists. It is seen that in the same sense as the sharp fc-space filter seems to be an idealistic choice to solve the 
mass function problem, the related fundamental stochastic diffusion process with its vanishing drift term must also be 
regarded as idealistic. 

No consistent solutions are, however, attainable for A\ ^ in combination with the original choice in the excursion 
set theory for the quasi time variable A = a 2 , assuming that a 2 is defined via a process-specific filter function and 
the power spectrum in the standard way. The inconsistency is seen by choosing, for example, E{A/|/} = Af in Eqs. 
(3.4) of Bond et al. (1991), where A is a constant independent of / and A, that is, Ax — A f in ^j. For A 2 = D, this 
is the order-1 or Ornstein-Uhlenbeck stochastic diffusion process. Appendix A summarizes some basic properties of 
this diffusion process. Assuming A > 0, as A — > 00 the variances (see Eq. A. 5 ) tend to the finite value, a 2 = D/(2A), 
but in the excursion set theory it is assumed that a 2 = A which tends to infinity. This contradiction shows that the 
specific choice, A = a 2 , is not suitable for a discussion of more general stochastic diffusion processes. 

However, a general resolution scale relation, A(R), is needed to handle different diffusion processes in an homoge- 
neous way. The most natural choice is the A(R) scale relation defined by the fundamental process, 
j_ 

A (^) = A l R dkk 2 P(k) - -L R-(n+3) for p {k) = p Qk n ( „ > _g) (1Q) 

as the common resolution scale for all diffusion processes investigated. A further motivation for this specific choice 
is given by the transformation property of the Wiener process. For example, the Wiener process, df(A) = dw(A), 
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is converted into the Ornstein-Uhlenbeck process, df(A) = —AfdA + \/Ddw(A), by the transformation /(A) = 
h(A)f w [T(A)}, with the amplitude transformation h(A) — e~ AA , and the resolution variable transformation T(A) = 
D I (2A)e 2AA , where again - and this is the important point - the resolution variable of the fundamental process, A, sets 
the basic resolution scale. Notice that the instrumental distinction between A and a 2 is necessary because A ^ a 2 in 
general. Only for the fundamental (Wiener) process we have A = a 2 . In the following sections pseudo mass functions, 
g(A, / c , /o), are derived and transformed into 'true' differential mass functions, n(M), using the Eqs. (|l^, |l6| ). 

6. Mass functions from order-0 (Wiener) stochastic diffusion processes 

Order-0 Gaussian stochastic diffusion processes are characterized by the constants A\ = —fx (or A = —fj, and v = in 
Eq. ||) , and A2 — D in (j|) , used to write the Ito stochastic differential equation in the form 

df(A) = fidA + y/Ddw(A), (20) 
with the solution 

,A 

/(A) = /o + /iA + VD / dw(A') , (21) 







and hence the instantaneous mean E[/(A)] = /o + //A, and the instantaneous variance var[/(A)] = DA. The 
corresponding Kolmogorov's forward equation is 

d d 2 u du an , - , , 

2W~^ = 0A ("~</</c)- W 
The filtered density contrasts are restricted to values below the zero redshift critical overdensity, e.g., f c = 
(3/20)(127r)( 2 / 3 )fi[J 0185 for the spherical symmetric cosmological infall model (Gunn & Gott 1972, Bertschinger 1985, 
Navarro, Frenk & White 1996). Notice that in (|22|), II is the probability density for / conditional to /(0) = /o, and 
that the process sample paths do not reach f c in the resolution interval (0, A). The cumulative distribution function, 
F(fo, f c , A) = 1 — F(fo, fc. A), defined by the probability density of first absorbed sample paths, is directly related to 
the pseudo mass function via eqs. jl^), (|l3|), and jl5|). We thus have 

F(/o,/ c ,A) = 1 - / C U(f J > A)df, (23) 



and compute F(fo,f c , A) solving ( p2| ) under the boundary conditions 

n(/, A = 0) = Su(f) , n(/ = /c, A) = (A > 0) . ' (24) 
Whereas the first equation reflects the deterministic initial state of the diffusion process, the second equation gives the 
appropriate boundary condition for an absorbing barrier at f c (e.g., Cox & Miller, 1965, p. 219): Assume a nonvanishing 
infinitesimal variance when / = f c . The probability that the trajectory is absorbed within AA can be approximated 
by the product of the probability that the trajectory is near (i.e., slightly below) / c , and the probability that the next 
increment A/ - which is proportional to \J AA (see Sect.[|) - carries it beyond f c . However, the latter probability is 
assumed to be always larger zero. In order to keep the probability density (i.e., the absorption probability per AA) 
finite that the trajectory is absorbed when / —> / c , the former probability must be zero because vAA/AA — > 00 when 
AA^O. 
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Fig. 4. Pseudo mass functions from order-0 processes for different fi values (continuous lines). For comparison, the Press- 
Schechter pseudo mass function = 0) is given (dashed lines). Both the mirror image and the Laplace transform method give 
the same results. 



6.1. Pseudo mass functions 
The unconstrained case of Eq 

ni(/ ,/,A) 1 



V^kDA 



exp 



(|22| ) has the general solution (Fig. 

(/-/o-M) 2 " 

2DA 



(25) 



where the sample paths start at fo. Exploiting the symmetry of the diffusion process, which still holds for /j, ^ 0, 
a solution is obtained by placing one source of sample paths at fo = 0, and a second source at fo = 2/ c (mirror 
image method), that is, n(/,A) = IIi(0,/,A) + a ITi(2/ c , /, A), where the factor a is determined by the condition 
Il(/ C , A) = 0. The solution gives the probability of nonabsorbed sample paths and takes the form 



n(/,A) 



i 



V2ttL>A 





(/-M) 2 " 


jexp 


2DA 



exp 



D 



exp 



jf - 2/ c - M) 2 
2L>A 



(-00 < / < / c ) . 



(26) 



It is seen that the second source in (|26|), i.e., the last term on the right-hand side, has a negative amplitude and acts 
as a 'sink' for sample paths. The drift parameter, /i, introduces a mass (A)-dependcnt shift of III which increases 
(fi > 0) or decreases (/j, < 0) the probability for absorbed sample paths (see shaded regions in Fig. |). Although the 
threshold, f c , is assumed to be mass-independent, the systematic drift might be interpreted as an effective change of 
the critical threshold in a mass-dependent way. Notice that a possible physical interpretation is that for high mass 
(low-A) clumps the effect of the surrounding volume elements is small so that the collaps will be spherical along all 
three directions with f c close to the spherical collaps model (Bernadeau 1994). For small mass (high A) clumps the 
effects of the surrounding volume elements might be more significant and can thus alter the effective value of f c . 
Substituting ( p6| ) into (|l^) gives the pseudo mass function for the order-0 diffusion process, 



<?(A,/ C 



HfcA) 
dA 



fc 



exp 



2L>A 



(27) 



The mirror image method is widely used in solving problems of heat conduction and heat diffusion (e.g., Landau & 
Lifschitz 1966). It 'reduces' the solution of the diffusion problem originally restricted to a specific domain of diffusion 
space to an infinite region. 

Since the absorbing barrier is A-independent, exact pseudo mass functions can also be deduced by taking Laplace 
transforms in the Kolmogorov's backward equation (Laplace transform method) giving the partial differential equation 
for 3*(A,/c,/o), 

f i<f(A,/c,/o) + ^7^5*(A,/c/o) = Xg*(XJ c Jo). (28) 

The backward equation is the appropriate equation because of the fixed final state f c , we seek for solutions g* as a 
function of the backward variable fo and set /o = later. For order-0 diffusion processes the general solution of (28) 
is 

<f(A,/c,/o) = Aef° H ^ + BeA* a « . _ (29) 

Inserting ( p9| ) in ( |28| ) gives a quadratic equation for the £Ts, and thus, after applying the condition g*(\,f c ,fo) < 1 
and g* (X, f c , fa = fc) — 1 (immediate absorption) from which the constants A and B are determined, the Laplace 
transform of the pseudo mass function (for f = 0) 

'fc 

D 



g*(X,f c ) = exp 



(/i - vV + 2D A) 



(30) 
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Fig. 5. Filter functions in configuration space: the sharp fc-space filter (dashed lines), and the truncated exponential filter 
(continuous lines) for the filter scale, R = 0.25 Mpc. The drift parameters have values between zero (dashed line) and 0.3. The 
amplitude of the power spectrum is Pq = 200 Mpc, and the spectral index n = — 2 (normalization and shape similar to a CDM 
power spectrum). The normalized Hubble constant, h — 0.5, is used in units of Ho = 100 km s -1 Mpc -1 . 



Inversion of (J3J 
For fj, — > 0, Eq. 



again gives (^): both the mirror image and the Laplace transform method lead to the same result. 
|27|) is the well-known Press-Schechter pseudo mass function with the correct normalization (see Fig. Eh. 



6.2. Mass functions 

Equating the process variance, DA, of the /'s with the mass variance, <r 2 (R), and using ( |l9|) gives the relation which 

determines the filter function for this process, 

1 r°° D f'T> 

dkk 2 \W R {k)\ 2 P(k) = DA = — dkk 2 P{k). (31) 



2tt 2 J 2tt 2 ... 

As expected, the filter function is the sharp fc-space filter (Eq. 0) . The normalization of the filter function yields D = 1 . 
Note that the relation between A and a 2 is simply given by <?(A) = A. In configuration space, the volume and the 
mass within the filter radius R as well as the sharp /c-space filter function (Fig.|^) are 



V s = 



W R (r = 0) 
1 



= 6tt 2 R 3 . 



For a scale-invariant power spectrum, P(k) 



Ppfc™, the M(A) relation becomes 



(32) 
(33) 

2vr 2 (n + 3)Aj ' (34) 
The condition for hierarchical growth of structure (variance decreases with mass) is n > —3. Notice that A is basically 
determined by the filter scale R, and by the power spectrum of the mass density distribution. This has important 



67T 2 po 



Pn 



implications when mass functions, n(M), are derived by using Eq. (15). The A(M) relation (we omit the subscript of 
M), 

Po f^Po\^ (35) 



(36) 



A(M) = , . 
where p = 2.7755 x 10 11 fl h 2 M Q Mpc -3 , gives the mass function, 



n(M) 



187T 2 V2^ 



(n 



6 + n f r 

3)— — ££ 



M P " +3 



3-ti 

. A 6 + 2 " 



exp 



(/c - M) 2 



2A 



Ratios of the mass functions induced by order-0 processes with positive drift parameters to the Press-Schcchtcr mass 
function are shown in Fig.|^. Eq. (|36| ) differs from the standard Press-Schechter function by the extra pA term in the 
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Fig. 6. Ratios of the mass functions induced by order-0 (left) and order-1 processes (right). A scale-invariant power spectrum 
with the index n = —2 and the amplitude Pq = 200 Mpc, the critical density contrast, f c — 1.68, the normalized Hubble constant, 
h — 0.5, and the density parameter, Qo — 1-0, are used. The different curves are labeled by the value of the corresponding drift 
parameter. Masses are given in solar units. 

exponential. As expected positive drift parameters increase the fraction of high / and hence the number of high-mass 
objects to the expense of objects with lower mass. For // — > 0, Eq. (|36| ) converges to the Press-Schechter mass function. 



7. Mass functions from order-1 (Ornstein-Uhlenbeck) stochastic diffusion processes 

Order-1 Gaussian stochastic diffusion processes are characterized by the parameters A± — Af (or v = 1 in Eq. ||), and 
Ai = D, used to write the Ito stochastic differential equation in the form 

d/(A) = —Af dA + VDdw(A). (37) 
Historical notes on this equation can be found in Wax (1954). Some basic informations about this process are given 
in Appendix A. The formal solution for /o = is 

r A 

/(A) = y/D / e - A(A ~ A > dw{A') , (38) 

"'0 

yielding the instantaneous mean, E[/(A)] = 0, and the instantaneous variance, 

var[/(A)]=^(l-e- 2 ^). (39) 

It is straightforward to write down the Kolmogorov's forward and backward equations (|6||7|) for this process. The 
boundary conditions are identical to those already used for the order-0 processes in Sect. 



7.1. Pseudo mass functions 



Since the absorbing barrier at f c is constant, the Laplace transform method (Sect. 6.1) can be applied. The Laplace 
transform of the corresponding Kolmogorov's backward equation (e.g., Beekman 1975 and the references given therein) 
leads to the pseudo mass function (definitions of the specific mathematical functions used in the following can be found 
in, e.g., Abramowitz & Stegun 1979), 

AS 2 \ D_ x/A (0) $ A a(A,0) 



g*(X,fc) = exp - . 

2D J D_ x/A (-U^) ^(A,/c) 

where D a (z) is the parabolic cylinder function with a = —A/A 
1 « (2/c^) ^ fn + X/A 



4> 



D,A 



(A,/ c ) = 



r( A )r(i± AM, ^ 



■ r 



-) 



(40) 



(41) 



n! ' V 2 

2A> \ 2 > ra=0 v 

an auxiliary function, and F(-) the Euler gamma function. The inversion of the Laplace transform appeared recently 
in the mathematical literature (see Ricciardi & Sato 1988) and involves the zeros A = X-wifc) of ( f4l| ) - they are all 
negative - giving the pseudo mass function 



fl(A,/ c 



OO 



^ K p (f c ) exp[A p (/ c )A] 

p=0 



(42) 
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Fig. 7. Pseudo mass functions for order- 1 processes obtained with the Laplace transform method (continuous lines) and with 
the mirror image method (dashed lines). Left (A = 0.1): f c = 0.3 and 1.68. Right (/ c = 1.68): A = 0.10, 0.01, 0.0025. Note the 
different scales of the two figures. 



where 

71 = 



,1 



2.4 



n + X p /A 



- * 



1 + X P /A 



2A 



(43) 



and W(-) being the digamma function. The asymptotic behaviour of 
' ^(i) 1 e-/o/(2A) . A small 



|) may be summarized as 



g(A,f c ) ~ < vSfl u ^ ' . (44) 

9{ ' \ K (f c )e-\^ : A large 1 ' 

The pseudo mass function can thus be expressed as the sum of exponential functions with negative exponents (A p < 0) 

and reduces to single exponentials for small and for large A. For A — ► 0, Eq. (^) tends to the Press-Schechter pseudo 

mass function. 



T hese results are compared to the closed-form analytic approximation obtained with the mirror image method (see 
Sect. |6.l[). The unconstrained Fokker-Planck equation has the solution 



ni(/ ,/,A) 



i 



/2vr 



2.4 



D{1 



= -2AA 



exp 



Ajf - f e~ AA ) 2 
D(l-e- 2AA ) 



(45) 



where all sample paths start at /o . The superposition of two sources of sample paths gives the constrained solution 



n(/,A) 



i 



/2vr 



2A 



D(l 



= -2AA 



A(J-2f c 



D(l- 



Taking into account the boundary conditions ( |24|) we obtain 



5(A, / c 



fc 



2A 



3 -2AA 



exp 



Aft 



D{l-e- 2AK ) 



AA 



(46) 



(47) 



As for order-0 processes, A is given by (|T|). Note that for fixed / c all pseudo mass functions agree in the low- A range, 
i.e., for large masses; differences are seen in the high- A range, i.e., for low masses. For small f c or small A the solutions 
obtained with the Laplace transform method can be approximated remarkebly well by the closed-form analytic results 
obtained with the mirror image method although the strict symmetry of the diffusion process is distorted by the 
/-dependence of the drift coefficient. For small A both solutions tend to the Press-Schechter pseudo mass function 
(see Figs. and ||). 



7.2. Mass functions 

Equating the process variance ( J38| ) with the mass variance, o~ 2 (R), and using (|l^) gives the relation which determines 
the window function, 



2tt 2 



dk k 2 \W R {k)\ 2 P{k) 




exp 



.4 



dkk 2 P{k) 



(48) 



For a given diffusion process the filter function and thus the region where material ultimately collapses to form a 
virialized halo cannot be chosen in any case independently from the power spectrum of the Gaussian random field. 
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For scale-invariant power spectra the filter function consistent with the density field and diffusion process is 

A P a fc" + 3 S 



W R (k) 



exp 



( AP k" + 3 \ 
\ 2^(n+3) J 







k < l/R 
else 



(49) 



This is a particular solution of the singular integral equation (48) found 'by eye'. Possible ambiguous solutions in the 
form of strongly fluctuating window functions in fc-space are physically not very interesting and can thus be disgarded. 
The normalization, D = 1, ensures Wa(k = 0) = 1. For A > it is an exponential filter truncated at the wavenumber 
1 / R. For large R Eq. ( f49| ) resembles the sharp fc-space filter whereas for small R the wings of the exponential get more 
important. Consider first a power spectrum with n — — 2 (the more general case n > — 3 is treated thereafter). For 
this realistic case the truncated exponential filter (M9) transformed into configuration space is 



W R {r) 



2ir 2 r 



■AP 
. 2tt 2 



r 
R 



[2n*R (^f + r 2 
AP r 



r) 



r 2 ((^) 2 +r 2 ) 



cos 



AP r 62 



9-7T- n 



Eq. (^) corresponds to 
by the inverse of W R (r) 



@ as A^O (Fig.|). 
as r — > 0, resulting in 



(50) 



The cumulative mass and volume at the filter scale R arc determined 



M E = po Ve , 
with the limits 

Ve = | ^ 



AP 
2ir 2 



e2» 



4tt 2 



APn 



( Ap a \ 2 

^-1) -- 



AP 



2vr 2 i? 2 



6vr 2 R 3 
(AP ) 3 



R « APo 



The M(A) relation thus becomes 



M E = po 



2tt 2 



Poe 



AA 



2 

A 



1 -1 



- 1) - 2A - AA 2 



(51) 



(52) 



(53) 



which can be solved numerically for A to compute the mass function (Eq. |15| ) . Ratios of the mass functions induced 
by order-1 processes with positive drift parameters to the Press-Schechter mass function are shown in Fig.^|. Eq. (53) 
can be approximated by 



Mi 



4tt 4 V A 3 



6 J 



leading to the A(M) relation (we omit the subscript of M) 
/4vr 4 Af A 3 \~^ 

~ [mi* " TJ ; 

and with the approximation (p 



to the mass function 



n{M) 



4tt 4 



fc 



9^2^ MP 3 



A 4 



2A 



-2AA 



exp 



Afc 



,-2AA 



AA 



(55) 



(56) 



It is seen that the mass function gives more objects with low mass and similar numbers of objects with high mass 
compared to the Press-Schechter case (see Fig. |J). The differences increase with increasing drift parameter and illustrate 
the effects of the different shapes of the filter functions (see Fig.||) on the mass function. 
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For scale-invariant power spectra with index n > —3, expansion of (^9|) and pice-wise integration gives the general 
M (A) relation, 



M E (A) = p V E , V e 1 = W R (r = 0) = 



1 

2^ 



dkk 2 W R {k) 



m=0 



(-l) m A™ . r . 
m! [m(^) + lj 



67T 2 



27r 2 (n+3) 



(57) 



from which the mass function can be obtained numerically in the same way as from (|53|). For A — ► only the m = 
term in the numerator of (|57|) gives a contribution to the sum, thus reproducing the volume of the sharp k space filter. 



8. Summary and discussion 

The present paper concentrates on the development of methods to derive analytic mass functions for given stochastic 
diffusion processes. For practical purposes the analysis is restricted to order-0 and order-1 processes. Whereas order-0 
processes are characterized by drift coefficients which are independent of the filtered density contrast, order-1 processes 
have a linear dependence on the filtered density contrast. It is seen that in order-1 stochastic diffusion processes, 
although derived within the framework of Markovian processes, realistic nonzero covariances appear in a very natural 
way without sampling over the diffusion (merging) history. The correlation between two mass scales gets important 
especially for small masses. For realistic cases the filter functions corresponding to the diffusion processes are derived. 
Whereas for all order-0 processes the filter function is the sharp fc-space filter, for order-1 processes exponential filters 
are found. The effects of non-sharp fc-space filters on the mass function are in the low-mass range where an excess 
compared to the standard Press-Schechter expectation is found. Future investigations should work with higher order 
diffusion processes in order to get more realistic potential wells (determined by the drift term) for sample diffusion. Of 
special interest are asymmetric potentials which might give a more appropriate handling of the asymmetric distribution 
of / within [— l,oo). 

Nonetheless, do we really need a discussion of more general diffusion processes? From the practical point of view 
one might argue that the Press-Schechter mass function, although derived under very specific assumptions, is well 
confirmed by comparison with N-body simulations (e.g., Efstathiou et al. 1988, White, Efstathiou & Frenk 1993, 
Lacey &: Cole 1994, Mo et al. 1996) and should thus suffice to be used as a simple fitting formula. The mass resolution 
and dynamic range of the simulations are, however, too low to test the prescription in detail. New results show 
indeed significant deviations from the analytical expectation (e.g., Jing 1998, Sheth & Tormen 1999, Governato et al. 
1999). In the light of the discussion given above one might expect such deviations from the Press-Schechter formula. 
Consequently one should avoid any unnecessary conditions on the type of diffusion process. The new coefficients offer 
a more flexible summary of the new mass functions. 

In forthcoming papers we want to derive consistent mass and filter functions for more complicated stochastic 
diffusion processes and compare the results with detailed N-body simulations and observed mass and luminosity 
functions derived from physically well-defined X-ray selected samples of clusters of galaxies. 
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for many useful comments. PS. thanks for support by the Verbundforschung under the grant No. 50OR9708 35. 



Appendix A: Order-1 stochastic diffusion processes 



The Ito stochastic differential equation of these processes is 
d/(A) = -AfdA + -/Ddw(A). 



(A.l) 



Early versions of this equation can be found in Uhlenbeck & Ornstein (1930) and Doob (1942). Eq. (A.l) can be 
solved by introducing the variable u = fe AA , and using the Ito formula for a sufficiently smooth function, /, of /(A) 
satisfying (g), 



#[/»] 



dA 



M (/(A), A) 



df 



1 d 2 P 

j^(/(A),A)^ 



dA + 



A 2 (f(A),A)^-dw(A). 
df 



(A.2) 



where the partial derivatives are evaluated at (A, /(A)). For / linear in / we have d 2 f/df 2 = 0, and (A.2) reduces 
to the usual chain rule. If / = u this gives du = \f~D e AA dw(A), which can be integrated and yields for the initial 
condition, /(0) = 0, 



/(A) 



,-A(A-A') 



dw(A') . 



(A.3) 
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The increments, dw(A), of the fundamental (Wiener) process have a Gaussian distribution (Sect.||). Hence, the linear 
superposition in (|A.3|) shows that the / are again Gaussian distributed with zero mean, the covariance 



E[/(Ai) /(A 2 ) ] = D-E 



Ai 



,-A(A!-Ai) 



du7(A;; 



,-A(A 2 -A 2 ) 



dw(A' 2 ) 



and the instantaneous variance 



D 
2A 



,-^|Ai-A 2 | 



3 -A(Ai+A 2 ) 



/(A) 



D 
2A 



1 



-2AA\ 



,(A.4) 



(A.5) 



Nonzero covariances between /(A) and the increments A/(A + /A) lead to more realistic correlations between different 
parts of each sample path, and thus to more realistic shapes of the corresponding process filters. For A > 0, and 



Ai, A 2 » 1, the last term in (A. 4) can be neglected, and the process has, as expected, a stationary covariance function 
which decreases exponentially with increasing distance |Ai — A2I, corresponding to a non- white power spectrum. For 
A < less fast decreasing covariances are found. 

The process can be gen eralize d to the case with A-dependent coefficients, 
d/(A) = -A(A) / dA + y/D(A) dw(A) . ' (A.6) 

For /(A = 0) = the solution of (O) is 



/(A) 



dw(A') y/D{A') exp 



A(s)di 



giving the covariance 

/>min(Ai ,A9) 

E[/(Ai) /(A a )] = / dA'D(A') exp 



and the instantaneous variance 

var[/(A)] = f dA' D(A') exp 
Jo 



Ai 



A(s) ds 



A(s') ds' 



2 A(s)ds 



(A.7) 
(A.8) 
(A.9) 
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